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(^ ' function for "gravitational entropy." We calculate its 

^^ ■ behavior for gravitational and density perturbations in 

closed, open and flat cosmological models and find that 
in all cases it increases monotonically. We are also able 
to calculate the entropy density of gravitational radia- 
{3JQ, tion produced by inflation. We compare the results with 

the behavior of the Weyl-tensor squared. Applying the 
formalism to black holes has proved more problematical. 



It is appropriate that I speak today on the subject of gravita- 
tional entropy because it was George Ellis, on my last visit to 
Cape Town, who instigated my investigations into this area — 
though he shouldn't be held responsible for any of the results. 

The problem is well-known. Ordinary thermodynamic systems, 
a box of gas for example, tend to grow more homogeneous with 



time, whereas gravitating systems tend to become more inhomo- 
geneous with time. In this sense gravitating systems are "an- 
t it her mo dynamic." I'm sure you've all seen the picture in Roger 
Penrose's book that shows this behavior. Indeed, the tendency 
toward inhomogeneity of a gravitating system can be viewed as a 
manifestation of the long-range nature of the gravitational force, 
which tends to cause the components of the system to clump. 

Now, ordinarily, for thermodynamic system, we associate the 
increase of homogeneity with an increase in thermodynamic en- 
tropy. This is a sign convention. One could, after all, choose 
Boltzmann's original H as the entropy function, and it would de- 
crease with homogeneity. Whichever sign one chooses, though, 
"gravitational entropy" or the "gravitational arrow of time" 
points in the opposite direction to the thermodynamic arrow. 

The question then becomes: Can you find a generally covari- 
ant function that characterizes the tendency of the gravitational 
field to become more inhomogeneous with time, a "gravitational 
entropy" ? 

However, in tackling this problem one immediately encounters 
a major conceptual difficulty: general relativity is a dynami- 
cal, not a thermodynamical theory, and it does not deal with 
temperatures. As prima facie evidence, I have been keeping a 
rough count of how many times the word "temperature" has 
arisen in the past four days and I would say five, ±1, and in all 
cases it has referred to the cosmic microwave background. An- 
other way of stating the difficulty is that general relativity deals 
with far-from-equilibrium systems for which a temperature is 
not well-defined. One usually, for example, doesn't talk about 
the temperature of a pendulum. Of course, you could define an 
effective temperature l/2mv'^ ~ kT but the fact remains- 

Jurgen Ehlers: The pendulum is not in equilibrium. 



TR: Exactly. 

So the question remains whether the concept of entropy can 
be incorporated into GR and how to do it. There don't seem 
to have been too many attempts in the hterature to define a 
gravitational entropy. The most well-known is Roger Penrose's 
Weyl-tensor hypothesis [0, the idea being that the square of 
the Weyl-tensor should be zero at the Big Bang but increase 
monotonically thereafter. However, the work of Bonnor [Q], as 
well as that of Wainright and collaborators [^ ^, ||, have thrown 
some doubt on the proposal. It's not my intent to criticize the 
hypothesis but rather to rethink the problem. Smolin [§ and Hu 
and Kandrup ||^ have examined some aspects of gravitational 
entropy, but not the arrow-of-time question. After I began work 
on this topic a paper by Brandenburger et al. appeared [Q, 
which again is concerned with the thermodynamical aspects of 
the entropy of gravitational radiation, but in the places where 
we can compare our work with their's, it seems to agree. Most 
of this talk will be based on a paper I wrote with Peter Anninos 
that appeared in Phys Rev D 0. 

Roger Penrose: By way of clarification I should say that I never 
meant the Weyl tensor to be a measure of gravitational entropy. 
I merely wanted it to be zero at the big bang. 

TR: Really? Well, if that's true I apologize. 

Roger Penrose: You're not the only person to have that misun- 
derstanding. 

TR: It does seem to be widespread. Perhaps we can sort out 
where people got this idea. 

At any rate, I decided to tackle the problem by a direct statistical- 
mechanics approach. We examine the definition 

S = \nQ (1) 



where S = gravitational entropy and Q = volume of phase space 
for the field. 

Why choose such a definition? In thermodynamic systems the 
phase space is much too hard to calculate. Instead you generally 
evaluate the partition function Z = Y^ie~^^\ from which the 
entropy is readily derived(here /3 = 1/(T) and E = the mean 
energy.) But this requires a temperature, which we don't have. 

On the other hand, consider a collection of A^ harmonic oscilla- 
tors with the usual Hamiltonian 

1 TV > TV 

H = ^T.'P^+^T.^' (2) 

This has a phase space. In fact, for one pendulum you know the 
phase space: an ellipse. If you increased the energy of the pendu- 
lum it would describe a slightly larger ellipse and the logarithm 
of the area between the two ellipses could be taken as entropy. 
However,we don't talk about entropy of dynamical systems-you 
know exactly the trajectory of the pendulum and so it doesn't 
make sense to define the logarithm of the area as entropy. You 
need an ergodic system, some lack of information. 

However, if we consider a system of oscillators with random 
phases, then we don't know where the oscillators are in this 
2N-dimensional phase space and it apparently makes sense to 
interpret the logarithm of this volume as entropy. In essence 
I am going to model the field as a collection of oscillators and 
treat the oscillators as a microcanonical ensemble, in which you 
assume any region of phase space is occupied with equal proba- 
bility. In that case S" = Infi is equivalent to S = —J2i Pi Inpj. 

To give you an idea of how it works, for the Hamltonian given 
in (0) above, we can actually calculate the phase-space volume 
analytically, using beta-functions, much in the way you'd calcu- 
late the volume of a sphere in elementary calculus. The result 



a - P''"g" (3) 

Noting that H = R^, this indeed is the volume of a 2N-dimensional 
sphere, or more precisely an ellipsoid. Taking S = InVt ^ A^ one 
can recover from this formula in the appropriate limits : 

• Entropy of an ideal gas(V = A; = 0) 

• Classical limit of entropy of an Einstein solid {H = E = NT) 

This isn't too surprising, as Einstein modeled a solid as a col- 
lection of harmonic oscillators. You can also get approximately 

• Entropy of the electromagnetic field. 
In fact you can write down 

• Partition function (|^) 

although I haven't figured out what to do with it. 

You might well be asking what this has to do with the gravita- 
tional field. The first reason to examine the above definition of 
entropy is that it does not require a temperature. The second 
reason is that there exists a Hamiltonian formulation of gravity 
(the ADM formulation). One can write down Hamiltonians for 
gravitational systems and thus one should be able to define Q 
below H. 

Here one needs to make an important caveat: In general, H 
will be time dependent, and therefore not the energy. So this 
is a slight extension of the usual statistical mechanics definition 
of energy. This doesn't necessarily bother me. In stat mech 
you for a system in contact with a heat reservoir the change in 
entropy is A^ = AQ/T. If you believe that Q is the kinetic 
energy of the Hamiltonian of the system, this does assume a 
change in H, although a quasi-static one. A time-dependent H 



is characteristic of an open system, which is really what we have 
in GR, and it is this time-dependence that will cause a change 
in entropy. So how do we calculate the change in HI 

Well, finally, I'll write down a metric, one for gravitational-wave 
perturbations: 



ds^ = a?{r]) —dif + \5ij + hij{ri, z)]dx^dx^ , (4) 

Here rj = conformal time; a{ri) = expansion scale factor; 

and hij <^ 6ij = represent the gravitational wave perturbations. 

Assuming a single polarization, the Einstein-Hilbert action to 
second order in the perturbation variables is: 

^ ''a^{h^-h'^)d'x, (5) 



647r 
Here (•) = d/drj and (') = d/dz. 

It's sometimes useful, in particular when comparing with other 



work, to change variables to 
above Lagrangian density: 



ih/y/32n, and you get for the 
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a 
a 



(6) 



The corresponding Hamiltonian density is by definition 7i 
nq — C =, or in this case. 



n 



02 + 0/2 _ ^^2 

a 



(7) 



The Hamiltonian is H = J Tid^x. In our case we want to make 
contact with the system of N harmonic oscillators and so we 
discretize the integral: 

1 ^ n 

H=^T.{-1+^?-^^), (8) 



where ¥ = 7i{L/Nf/^, = ^{L/Nf/'^, H = H/N^, and L is an 
arbitrary length scale. 

Notice that this equation resembles the harmonic oscillator Hamil- 
tonian Eq @ except for two features. The first is that it contains 
a gradient term: 0' = 0j — 0(i+i) . This is not a big problem. By 
defining a new variable ^j = 0j — 0(i+i), one can integrate this 
"nearest-neighbor potential" in the same way one does the har- 
monic oscillator potential and get an result for the phase space 
identical to Eq (^), but for insignificant numerical factors. 

The second feature is the main problem: the minus sign in front 
of the harmonic-like potential of the last term. This implies, 
first, that H can change sign and that, second, we have a re- 
flection barrier instead of a potential well. This implies that the 
phase space is unbounded. Without going into the gory details, 
I claim that in the perturbative limit, one can still analytically 
integrate the phase space using hypergeometric functions (see 
Rothman and Anninos for specifics), and one formally recovers 
expression (|^). Nevertheless, one is forced to impose some sort 
of cutoff, in momentum space if // > or in configurations space 
ii H < 0. (It turns out that H > represents growing modes 
and H < represents decaying modes, which is an interesting 
result in its own right, but we do not go into details here. The 
growing modes, which result in an increase in entropy, are of the 
main concern.) There are several conceivable methods to impose 
the necessary cutoffs, but it turns out they all give qualitatively 
the same results. 

So, the procedure to calculate the entropy is to regard Q as 
constant on each spacelike hypersurface, in which case we can 
take d/a = k, the spring constant. To find Q, we first need H 
, which we get by considering the equations of motion resulting 



from the action (|^): 

h + 2-h- h" = 0. (9) 

a 

We assume a separable solution with random phases a^: 

Y.(h, + ^h, - /.;') = (h + ^h- h") E e-^ = 0, (10) 

j j 

The subscript j refers to different waves, not to different coor- 
dinates. By taking the phases to be random we are essentially 
assuming an incoherent source, and since the aj do not enter into 
the solutions, I supress them from now on. Eq (^ is Bessel's 
equation and for the matter-dominated epoch with a = aorf-, 
one gets: 

/iocr3/V±3/2(A;r/)e*'=^ (11) 

Now, the J's have standard asymptotic forms for A;?7 ^ 1 (A ^ 
Hubble radius )and krj^ 1 (A <^ Hubble radius). 

For krj <^ 1, we have: 



h 



h{kr])-' + h2e''', (12) 



where hi and /i2 are constants representing the decaying and 
growing modes respectively. (In this case the "growing" modes 
are constant in time.) 

For ki] ^ 1 we have: 

h oc Y ^ (^j cos{kr]) + sin{krj) e^^^ 

oc {kri)^"^ X [oscillations] (13) 

In terms of </> ~ ah, for kr] ^ 1: 

= Harmonic oscillator potential at fixed time, (14) 



and 

fl oc , oc constant x [oscillations]. (15) 

Therefore ^2 oscillates, but otherwise does not change. Given 
that in this limit if is a harmonic oscillator, this isn't too sur- 
prising. One could average over several cycles to obtain a strictly 
constant H: 

(^)7^= /rf^x[7r2 + A;V- VI- (16) 



In the nonlinear regime, Q would increase, which is encouraging 
for the interpretation of In Q as entropy. 

In the opposite limit kt] '^ 1, the gradient term in (|]) is negligi- 
ble and we get for the time dependence of if : if oc vr^ — acj)^ /a 
and n oc H^/{-d/a)^/\ 




andf^oc<' ''_ ' for growing modes, ,^^. 
1 ^ SAT ^ for decaying modes. 




This result is at first glance somewhat perplexing. Superhorizon 
modes should be frozen in, no oscillations allowed, so why should 
we get an increase in ^2? In terms of h one sees that indeed: 
(j) = ah + ah = ah. That is since h = 0, the superhorizon modes 
are frozen in and the increase in Q is due entirely to expansion 
d. Furthermore, since the superhorizon modes are frozen, we 
know the position of each oscillator and the assumption of ran- 
dom phases breaks down. This suggests that our definition of 
entropy is only apphcable to subhorizon processes, which is not 
necessarily a bad thing. Such an interpretation agrees with the 
conclusions of Brandenburger el al. 

Without going into further details, I merely state that one can 
repeat the above analysis for density perturbations (radiation 

9 



and dust) in flat, open and closed cosmologies. One always 
has to worry about gauge problems in inhomogeneous models 
and so we use the gauge-invariant formalism of Mukhanov et 
al. [|l^]. We find that in all cases Q nionotonically increases 



for growing modes (H > 0)and decreases for decaying modes 
(H < 0). This is of course encouraging for the interpretation of 
In fi as a measure of gravitational entropy. 



In one of George's preprints [11] he and his coauthor Reza Tevakol 
mentioned that it is hard to see how to apply this formalism to 
GR in full generality since the number of modes, N, is infinite 
in field theory. I should clear up this misunderstanding. We 
are not doing a second-quantized field theory. A^ is merely the 
number of Planckian oscillators- modes in a box-which should 
be finite if the box under consideration is of finite size. To be 
sure, N is merely the Jeans' number ~ u^a^, where u is the 
oscillator frequency. In all our models S ^^ N (see especially Eq 
(p6|) below), as is the case for electromagnetic radiation. 

From this fact, we can actually calculate the entropy density 
of gravitational radiation produced during inflation. I did this 
calculation quite recently and it is not to be found in Rothman 
and Anninos, but it does follow closely the method outhned by 
Peebles [Q on page 492-494. Our Lagrangian density can be 



rewritten slightly as 

C = rnlih^ - h'^/a") (18) 

where now ' = d/dt. This is a fairly generic Lagrangian density 
with m,ph playing the role of the usual kinetic energy cfp' term. 
As Peebles shows, mph ~ H, the Hubble constant, is the rms 
value of the field. If we assume that h represents the size of the 
perturbation during inflation, we then have 

'. ^ ^ ^ 5f 4. (IS) 
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where e is the energy scale of inflation. 

Now, as mentioned above, the entropy S ^ N. Define an en- 
tropy density a = N/a? = u^. But in our case, tu ~ if ~ niph 
and so 0" ~ h^rrip. All we need to do is scale this down to the 
present day to find cTq ~ hlrrip, where ho is today's value of h. 
Since u! ~ niph, we can also write ao ~ {niphoUJoY^'^- Now, hg is 
the strain that has just started oscillating as its corresponding 
proper wavelength Aq falls within the Hubble radius. Peebles 
shows that 

K - ny^XoHo^ (20) 



m; 



But Ao ~ l/ujo- Thus 



a, ~ n'J'H!/^4H (21) 



m 



3/2 



P 



With Ho ~ T^/ml, we get 



.3 



^o ~ n^T!^. (22) 



This says that the entropy of gravitational radiation produced 
during inflation is supressed from that of the microwave back- 
ground by the factor e^ /rrip. Apart from the factor f2^/^, which 
is of order unity, Eq (^2|) is identical to an equation in Branden- 
burger et al., obtained by much more difficult means. 

To conclude this part of the talk, we find that the candidate 
gravitational entropy function S* = Infi behaves reasonably on 
subhorizon scales in that it increases monotonically, and when 
applied to the case of gravitational radiation during infiation 
produces the result that a generic field theory should produce. 
Nonlinear calculations should be carried out. 
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I would like now to briefly discuss two other aspects of this 
project. We did compare the behavior of S with that of C^ for 
the metric 

ds'^ = a^{ri)[-{l + 2^{ri, z))dri'^ + {l-2<!?{ri, z))-fijdx'dx^], (23) 

where $ is the gauge-invariant version of h with solution $ = 
(mi +U2i]~^) e''^^^\ Specifically, we examined 

<-a/3 <-7'5 - ^^4 ' ^^^^ 

where the solution for $ is the one just given. I bear in mind 
Roger's earlier comment that he did not mean for C^ to be a 
measure of entropy, but with a r^ rf,we find that (|^) decreases 
with time and inhomogeneity, which seems a rather strange be- 
havior. You might think to correct this behavior by introducing 
an overall minus sign, but then the decaying modes increase with 
decreasing inhomogeneity, which is no less strange. 

On the other hand, we also examined the form of the Penrose 
hypothesis suggested by Wainwright and collaborators |P, |, |[ 
To lowest order this is: 



CjC^r _ 4a^ (<|.,..)^ _ r,^ ($,. 



i?„^i?°/5 9(d4 - aa^'d + a^'a^) 27 

which has a time dependence 

*-"a/3 *-"75 I Ti , tor growmg modes 



(25) 



oc 



Ra/sR"^ \ V 1 for decaying modes. 

Such a time dependence is not only reasonable but is,to one's 
great surprise (at least our great surprise), identical to that of 
the Hamiltonian H for k = dust. The expression (p5|) is an 
approximation, so we don't know whether we are facing a coin- 
cidence or something more profound. We haven't investigated 
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the matter further, but I think if there is a graduate student 
interested in numerical work, it might be worth looking into. 

I'd like now to turn to my recent-and so far not entirely successful- 
attempt to apply this formalism to black holes. Take expression 
(0) for the phase space of A^ harmonic oscillators. Now imagine 
creating a black hole out of quantum oscillators. In this case the 
Hamiltonian H should be the total energy,M and we should also 
have M = Ne, where e = average oscillator energy. Further the 
angular frequency uj = 2ne = k"^, where k is the spring constant. 
Plugging all this into (|^) and using Stirling's approximation on 
the N\ gives identically 

n = e^and S = N (26) 

Again , N = M/e = MX, where A is the wavelength. A priori, 
we expect A ~ AM. With this value we find 

S^4M^ = SBH/{27r^), (27) 

where Sbh is the canonical Bekenstein-Hawking value. 

About half the people I have shown this to call it a "discovery" 
and the other half call it a "dimensional coincidence." I think 
it is bit more than the latter. In his review article on string 



theory and black holes, Gary Horowitz |]14[ points out that string 
entropy is given by Sg ~ isMg, where ig is the length of the 
string and Mg is the string mass. In the limit that ig becomes 
the Schwarzschild radius and Mg becomes the mass of the black 
hole, we get S ~ RbhMbh- Horowitz calls this "remarkable," 
yet it seems to me to be essentially the same argument I just 
gave for ordinary oscillators. The result appears to be quite 
general. 

The main problem with this technique, as I see it, is that it 
doesn't give a value for A^, which must be put in by hand. The 
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question is, if you believe S" ~ A^, can you get a better fix on 
N7 The assumption that black-hole entropy is related to the 
number of interior modes has been challenged, but in light of 
the recent developments in string theory, it seems to me worth 
the effort to find a more modest description along these lines. 
What I've recently been attempting to do is count interior modes 
of a Schwarzschild black hole. 

The idea is fairly simple. Consider the usual scalar wave equa- 
tion on a Schwarzschild background: 

,^ 2M,^,^ 2M, , 2M , , 

- (1 ) (1 )^,rr + -^^A 

2M,,2M £(£+1), o, , , 

+{1 — + ^^ -^^^ =0 28 

This has a well-known effective potential (see MTW, p. 868 
||13|), which goes to minus infinity at the origin. But the only 
thing that happens aX R = 2M is that Veff goes to zero. More- 
over, although the Legendre-functions {P's and Q's)th.a.t are the 
solutions to (|28D in the limit of zero frequency blow up either at 
2M or infinity (the no-hair theorem), the P's are actually regu- 
lar on the interior. So the idea is to put in a Planck-frequency 
cutoff = Up and see if one can count modes in a sensible fashion. 



Now, you might think this is impossible, since ( pSD has no an- 
alytic solutions on the interior, but you don't actually need to 
solve the equation. Define a new variable related to the old by 
the integrating factor of the equation: ip = pAf-r -^-i/2^_ This 
eliminates the first derivative term in the wave equation. Re- 
calling that r is timelike for r < 2M, make the ansatz ip ~ e*'^'". 
Then you find 

f^^^.^m^f-Ml^^?E^/l±A]J,^o (29) 

Since t is spacelike for r < 2M, this is just like a flat-space wave 
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equation, 

i," + ki^;^, = 0, (30) 

where kl^j is the quantity in brackets in (^). 

Now you just count modes in a box. Following the suggestion 
of 't Hooft [|l^, the total number of radial modes beneath a 



frequency u will be, 



N^ = - keff didr (31) 

71 Jrl 



I claim that you can integrate this expression. The integrals 
converge aX R = 2M and the potential terms don't contribute. 
The only thing that contributes is the cutoff frequency and the 
radial cutoff (= 2M). In fact, it is quite remarkable that of the 
six terms in the final integral, the only non-negligible term is the 
one with the correct scaling; and the result is S' = (3/8)7ra;^M^. 
Choosing the Planck frequency Up = 2tt as the cutoff, the ex- 
pression becomes 5* = {3/16)7iSbh, or about .6 the canonical 
value. In flat space, one gets exactly .25 Sbh- 

Now, a scalar wave has only one polarization. If one counted 
gravitational perturbations instead, one would use the Zerilli or 
Regge- Wheeler equation, which have different potentials. But 
since the potential terms don't contribute, one gets the same 
answer per mode. Therefore, for the two modes of gravitational 
waves, the result is 1.2 Sbh, which seems too close for comfort. 

Unfortunately, there are two flies in the ointment. The first 
is that the integral over r in (|3l|) is a timelike integral, which 
makes the interpretation unclear. More seriously, I have ignored 
azimuthal modes. To count them, the integrand in (|3T|) should 
contain a factor {2i + 1). However, in that case the scaling goes 
wrong; i.e., one gets S ~ w^r^. 

I am at this time unable to justify ignoring the azimuthal modes 
and so I must let the result stand as what may well be a dimen- 
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sional fluke. Even so, it contains some interesting points in com- 



mon with earlier work of Bombelli et al. ||T6[. They calculate 
the entropy due to a system of coupled oscillators. To carry out 
the calculation they define a density matrix and perform partial 
traces over an imaginary sphere that separates "inside oscilla- 
tors" from "outside oscillators." In the continuum limit they 
find that the entropy of the oscillators within this imaginary 
sphere scales as the area. Although the considerations by which 
they obtain their results are somewhat different than mine, they 
nevertheless find that the result is independent of the potential 
and depends only on the cutoffs in uj and r and is proportional 
to the number of fields under consideration. (This will be true 
whether azimuthal modes are considered or not.) Whether this 
points to some deeper understanding of black hole entropy or is 
merely a reflection of the fact that both calculations are based 
on oscillator systems, or whether it is a basic reflection of the di- 
mensionful constants that govern the problem is something that 
I have yet to sort out. Thank you. 



References 

[1] Penrose, R. (1989), in Fourteenth Texas Symposium on 
Relativistic Astrophysics, ed. E. J. Fenyves (New York 
Academy of Sciences, New York). 

[2] Bonnor, W.B. (1987) Physics Letts. A 122, 305. 

[3] Wainwright, J. and Anderson, P.J. (1984), Gen. Rel. Grav. 
16, 609. 

[4] Goode, S.W. and Wainwright, J. (1985), Class. Quantum 
Grav. 2, 99. 



16 



[5] Goode, S.W., Coley, A.A. and Wainwright, J. (1992), Class. 
Quantum Grav. 9, 445. 

[6] Smolin, L. (1985) Gen. Rel. Grav. 17, 417. 

[7] Hu, B. and Kandrup, H. (1987) Phys. Rev. D 35, 1776. 

[8] Brandenberger, R., Mukhanov, V. and Prokopec, T. (1993) 
Phys. Rev. D 48, 2443. 

[9] Rothman, T. and Anninos, A. (1997) Phys. Rev. D 55, 
1948. 

[10] Mukhanov, V., Feldman, H. and Brandenberger, R. (1992) 
Phys. Rep. 215, 203. 

[11] Tavakol, R. and EUis, G.F.R., "On Gravitational Entropy," 
University of Cape Town, preprint, (1997). 

[12] Peebles, P.J.E., Principles of Physical 

Cosmoiog7(Princeton: Princeton University Press, 1993. 

[13] Misner, C.W., Thorne, K.S. and Wheeler, J.A. (1973), 
Gravitation, (W.H. Freeman and Company, New York). 

[14] Horowitz, Gary, |gr-qc/ 970407^ 

[15] 't Hooft, G. (1986), Nuc. Phys. B 256, 727. 

[16] BombeUi, L. et al. (1986), Phys. Rev D 34, 373. 



17 



